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by including quarks. 
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1. Introduction 

Past LGT simulations of the deconfining transition focused primarily on boundary conditions 
(BCs), which are favorable for reaching the infinite volume quantum continuum limit (thermo- 
dynamic limit of the textbooks) quickly On N t N 3 lattices these are periodic BCs in the spatial 
volume V = (aN s ) 3 , where a is the lattice spacing. The temperature of the system is given by 

' ^, (N T <N S ). (1.1) 



aN T L x 

In the following we set the physical scale by 

r c = 174MeV, (1.2) 

which is approximately the average from QCD estimates with two light flavor quarks. This implies 
for the temporal extension 

L T = aN T = 1.13 fermi. (1.3) 

For the deconfinement phase created in a heavy ion collision the infinite volume limit does not 
apply. Instead we have to take the finite volume continuum limit 

N S /N T = finite , N x — > °° , L T finite , (1.4) 

and periodic BCs are incorrect, because the outside is in the confined phase at low temperature. 
E.g., at the BNL RHIC one expects to create an ensemble of differently shaped and sized deconfined 
volumes. The largest volumes are those encountered in central collisions. A rough estimate of their 
size is 

71 x (0.6 x Au radius) 2 x c x (expansion time) 

= (55 fermi 2 ) x (a few fermi) (1.5) 

where c is the speed of light. Here we report on our work [jjj], which estimates such finite volume 
corrections for pure SU(3) and focuses on the continuum limit for 

L, = aN s = (5 - 10) fermi . (1.6) 

2. Equilibrium with Unconventional Boundary Conditions 

Statistical properties of a quantum system with Hamiltonian H in a continuum volume V, 
which is in equilibrium with a heatbath at physical temperature T, are determined by the partition 
function 

Z(T, V) = Tve- H ' T = £(0 \e- H ' T \<j>) , (2. 1) 

<P 

where the sum extends over all states and the Boltzmann constant is set to one. Imposing periodic 
boundary conditions in Euclidean time T and bounds of integration from to 1/r, one can rewrite 
the partition function in the path integral representation: 

Z{T,V)= J D<j>exp{- dTL E (<j>,<j>)\ . (2.2) 
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Nothing in this formulation requires to carry out the infinite volume limit. 

In the following we consider difficulties and effects encountered when one equilibrates a hot 
volume with cold boundaries by means of Monte Carlo (MC) simulations for which the updating 
process provides the equilibrium. We use the single plaquette Wilson action on a 4D hypercubic 
lattice. Numerical evidence suggests that SU(3) lattice gauge theory exhibits a weakly first-order 
deconfining phase transition at some coupling p t 8 (N t ) = 6/gj(N t ). The scaling behavior of the 
deconfining temperature is 

T c = c T A L (2.3) 

where the lambda lattice scale 

aA L = f k (p) = X(g 2 ) (bo8 2 )- bim) e- l/i2b0g2) , (2.4) 

has been determined in the literature. The coefficients bo and b\ are perturbatively determined by 
the renormalization group equation: 

11 3 34 / 3 \ 2 

' = Tl6^ and ' 1 = Tlj6^J • (2 " 5) 

Relying on work by the Bielefeld group [|2|] we parametrized in [J3j] higher perturbative and non- 
perturbative corrections by 

X(g 2 ) = \+a ie - ailgl + a 3 g 2 +a 4 g 4 with (2.6) 



a\ = 71553750, a 2 = 19.48099, a 3 = -0.03772473, a 4 = 0.5089052, which turns out to be in 
good agreement with [Q] in the coupling constant range for which the latter is claimed to be valid. 

Imagine an almost infinite space volume V = L%, which may have periodic BCs, and a smaller 
(very large, but small compared to V) sub-volume Vb = L% Q . The complement to Vb in V will be 
called V\ . The number of temporal lattice links N T is the same for both volumes. We denote the 
coupling by j8q for plaquettes in Vb and by j3f for plaquettes in Vi . For that purpose any plaquette 
touching a site in Vi is considered to be in Vf. This defines a BC, which we call disorder wall. 

We would like to find couplings so that scaling holds, while Vo is at temperature T = 174 MeV 
and Vi at room temperature 7\. Let us take j8f = 5.7 at the beginning of the SU(3) scaling region. 
We have 

i n 10 ~ Z° - 51 - j^Ml n i\ 

~Trao-fx(J>§) ( ° 

where a,- is the lattice spacing in Vj, i = 0, 1. Using the lambda scale yields j8q 25 and T c estimates 
of the literature give L T > 10 11 a. In practice we can only have j3 ( f in the scaling region. We keep 
up the relation 

^4 = - « 10 10 (2-8) 
(g/ai) a 

where £ is a correlation length. Therefore, %ja\ is very small and the strong coupling expansion 
implies j3f lO" 10 ' , i.e., j8f = 0. 
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Figure 1: Fits of pseudo-transition coupling constant values (left). Estimate of finite volume corrections to 
T c (right). 



3. Monte Carlo Calculations with Disorder Wall BCs 

In the disorder wall approximation of the cold exterior we can simply omit contributions from 
plaquettes, which involve links through the boundary. Due to the use of the strong coupling limit 
for the BCs, scaling of the results is not obvious. 

We use the maxima of the Polyakov loop susceptibility 

Zmax = ^[(| J P| 2 )-(| J P|) 2 ] max , J P = I^ (3-1) 

s x 

to define pseudo-transition couplings ^p t {N s ;N T ). For periodic BCs they have a finite size behavior 
of the form ^ 

P$,(N S ;NJ = P t s (N T )+a p 3 f^j + .... (3.2) 
Our BCs introduce an order iV? disturbance, so that 



The left Fig. [| shows thus obtained fits of pseudo-transition coupling constant values versus N z /N s 
(the N s — > °° value is extrapolated from simulations with periodic BCs). Using the scaling relation 
(2^) we eliminate the coupling in favor of T c and L s and obtain the right Fig. [[]. There are no 



free parameters in this step, because the scaling relation was determined previously in independent 
work. The N t = 4 and N t = 6 data collapse to one curve, i.e., despite the small values of the 
temporal lattice sizes the results are perfectly consistent with scaling. 

The left Fig. |2| shows the Polyakov loop susceptiblity on a 4 x 16 4 lattice with disorder BCs 
and its full width at 2/3 maximum, which we used instead of the more conventional full width at 
half maximum, because the former is easier to extract from MC data (smaller reweighting range). 
Our width data are fitted to the form 



2/3-1 \NJ 2 \N S 
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for periodic BCs and to 



^ =««(£)'+««(£] ™ 



for disorder wall BCs. The first term reflects in both cases the delta function singularity of a first 
order phase transition, i.e., the width times the Polyakov loop maximum is supposed to approach a 
constant for N s — ► °°. The leading order correction to that is 1/Volume for periodic BCs and 1 /N s 
for disorder wall BCs. Plots of the corresponding fits are shown in Figs. ^| (right) and || (left). 
As before, we use the scaling relation ( 2A ) to eliminate the coupling constant and show in Fig. [| 



(right) the thus obtained volume dependence of the width of the transition. Again, we see collapse 
to a nice scaling curve. 



4. Shortcomings of the Disorder Wall BCs 

The spatial lattice spacing a s should be the same on both sides of the boundary, but for the dis- 
order wall this is not true. It reflects the temperature jump at the price of introducing a similar jump 
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in the spatial lattice spacing. Its main advantage that it allows for technically simple simulations, 
and one can hope that the temperature jump is the only relevant quantity for the questions asked. 

A construction, called confinement wall in for which the physical length of one spacelike 
lattice spacing stays constant across the boundary can be achieved by using an anisotropic lattice 
for the volume V\ : 

S({U}) = f £ReTr(C/ a ) + ^£ReTr(C/ Df ) . (4.1) 

The lambda scale of this action has been investigated by Karsch [j|] and in the continuum limit one 
finds 

PI /Pf = (a s /a t ) 2 . (4.2) 

When we aim at «o = a s sa 10~ 10 a T the sublattice V\ is driven to /3f = and the simulation of 
the confined world becomes effectively 3D. However, in a first step one may be content with a 
temperature slightly below T c on the outside, so that the confinement wall allows to have all j8- 
values in their scaling regions. Another approach may want to rely on symmetric lattices to model 
low temperatures. 

5. Summary and Conclusions 

1. As noted before [Q] finite size corrections to deconfinement properties of SU(3) are very 
small for periodic BCs. 

2. For volumes of BNL RHIC size the magnitudes of SU(3) corrections due to cold boundaries 
appear to be comparable to those of including quarks into pure SU(3) LGT. Our data show 
the correct SU(3) scaling behavior. 

3. Extension of measurements should be done, to calculate the equation of state. 

4. Previous calculations ^ ^ of full QCD at finite temperatures and RHIC (low) densities 
should be extended to other than periodic BCs. 

5. There appears to be a variety of options to include cold boundaries and approaching the finite 
volume continuum limit. Therefore, more experience with pure SU(3) LGT is desirable be- 
fore including quarks. Next, we intend to focus on the confinement wall with both couplings 
in the scaling region (i.e., an outside temperature just below T c ). 
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